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Abstract. Low-mass dilepton spectra as measured in high-energy heavy -ion collisions are a unique 
tool to obtain spectroscopic information about the strongly interacting medium produced in these re- 
actions. Specifically, in-medium modifications of the vector spectral function, which is well known 
in the vacuum, can be deduced from the thermal radiation off the expanding QCD fireball. This, in 
particular, allows to investigate the fate of the p resonance in the dense medium, and possibly infer 
from it signatures of the (partial) restoration of chiral symmetry, which is spontaneously broken 
in the QCD vacuum. After briefly reviewing calculations of thermal dilepton emission rates from 
hot QCD matter, utilizing effective hadronic theory, lattice QCD or resummed perturbative QCD, 
we focus on applications to dilepton spectra at heavy-ion collider experiments at RHIC and LHC. 
Qh This includes invariant-mass spectra at full RHIC energy with transverse-momentum dependencies 

and azimuthal asymmetries, as well as a systematic investigation of the excitation function down to 
fixed-target energies, thus making contact to previous precision measurements at the SPS. Further- 
more, predictions for the energy frontier at the LHC are presented in both dielectron and dimuon 
channels. 
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The exploration of matter at extremes of temperature (T) and baryon density (pg) is 
at the forefront of research in contemporary nuclear physics, with intimate connections 
to high-energy, condensed-matter and even atomic physics [1]. Theoretical efforts over 
the last few decades are suggesting an extraordinary richness of the phase diagram of 
strongly interacting matter, which should ultimately emerge from the underlying theory 
of Quantum Chromodynamics (QCD) as part of the Standard Model. However, several 
basic questions, both qualitative and quantitative, such as the possible existence of first 
order transitions and their location as function of baryon- chemical potential (/j^) and 
temperature, remain open to date (see, e.g., Ref. [2]). A close interplay of experiment and 
theory is needed to create a robust knowledge about QCD phase structure. On the one 
hand, naturally occurring matter at temperatures close to or beyond the expected pseudo- 
critical one, r pc ~ 160 MeV [3, 4], may last have existed ~14 billion years ago, during 
the first tens of microseconds of the Universe. On the other hand, at small temperatures, 
matter with baryon densities close to or beyond the critical one for the transition into 
quark matter may prevail in the interior of compact stars today, but its verification and 
exploration from observational data is challenging [5]. It is quite fascinating that tiny 
man-made samples of hot QCD matter can nowadays be created and studied in the 



laboratory using ultrarelativistic heavy-ion collisions (URHICs). Significant progress 
has been made in understanding the properties of this medium through analyses of 
experiments conducted at the CERN's SPS, BNL's Relativistic Heavy-Ion Collider 
(RHIC) and CERN's Large Hadron Collider (LHC) (see, e.g., the recent Quark Matter 
conference proceedings [6, 7]). For example, systematic investigations of the produced 
hadron spectra have revealed a hydrodynamic behavior of the bulk matter in the region 
of low transverse momenta (q t < 2-3 GeV) and a strong absorption of hadron with high 
transverse momentum (q t >6GeV). Even hadrons containing a heavy quark (charm or 
bottom) exhibit substantial energy loss and collectivity due to their coupling to the 
expanding fireball. While the total charm and bottom yields are essentially conserved, 
the production of heavy quark-antiquark bound states (charmonia and bottomonia) is 
largely suppressed. The relation of the above hadronic observables to spectral properties 
of the medium is, however, rather indirect. Low-mass dileptons, on the other hand, are 
radiated from the interior of the medium throughout the fireball's lifetime, as their mean- 
free path is much larger than the size of the fireball. Thus, their invariant-mass spectra 
directly measure the in-medium vector spectral function, albeit in a superposition of the 
varying temperature in the fireball's expansion. 

The dilepton program at the SPS has produced remarkable results. The CERES/NA45 
dielectron data in Pb-Au collisions, and particularly the NA60 dimuon spectra in In-In 
collisions, have shown that the p -meson undergoes a strong broadening, even complete 
melting, of its resonance structure, with quantitative sensitivity to its spectral shape, see 
Refs. [8, 9, 10] for recent reviews. The QCD medium at SPS energies is characterized by 
a significant net-baryon content with chemical potentials of \Xb — 250 MeV at chemical 
freezeout, T c ^ ~ 160 MeV [11], and further increasing as the system cools down [12]. 
Baryons have been identified as a dominant contributor to the medium modifications of 
the p's spectral function. The question arises how these develop when moving toward 
the net bary on-free regime in the QCD phase diagram, /ig <C T . Theoretical expectations 
based on the hadronic many-body approach [13] suggest comparable medium effects in 
this regime, since the relevant quantity is the sum of baryon and anti-baryon densities, 
and this turns out to be similar at SPS and RHIC/LHC [12], at least close to T pc . 
Since T c ^ ~ T pc at collider energies, the total baryon density at RHIC and LHC in the 
subsequent hadronic evolution of the fireball will remain similar. We also note that the 
\Xb — OMeV regime is amenable to numerical lattice QCD calculations, both for the 
equation of state of the medium evolution, and in particular for the microscopic dilepton 
production rate, at least in the QGP phase for now [14, 15]. Furthermore, since the 
phase transition at /ig ~ MeV presumably is a continuous crossover, a realistic dilepton 
rate should vary smoothly when changing the temperature through T pc . Thus, after the 
successful fixed-target dilepton program at the CERN-SPS, the efforts and attention are 
now shifting to collider energies around experiments at RHIC and LHC. 

In the present article we will focus on the theory and phenomenology of dilepton 
production at collider energies (for a recent overview including an assessment of SPS 
data, see, e.g., Ref. [16]). The presented material is partly of review nature, but also 
contains thus far unpublished results, e.g., updates in the use of nonperturbative QGP 
dilepton rates and equation of state, and detailed predictions for invariant-mass and 
momentum spectra for ongoing and upcoming experiments at RHIC and LHC, including 



an excitation function of the beam energy scan program at RHIC. 

This article is organized as follows. In Sec. 2 we briefly review the calculation of 
the thermal dilepton emission rates from hadronic matter and the quark-gluon plasma 
(QGP). We elaborate on how recent lattice-QCD results at vanishing three-momentum 
(q=0) can be extended to finite q to enable their application to URHICs. In Sec. 3 we 
discuss in some detail the calculations of dilepton spectra suitable for comparison with 
experiment; this involves a brief discussion of the medium evolution in URHICs (includ- 
ing an update of the equation of state) in Sec. 3.1, and of non-thermal sources (primordial 
production and final-state decays) in Sec. 3.2. It will be followed by analyses of mass 
and momentum spectra, as well as elliptic flow at full RHIC energy in Sec. 3.3, and of an 
excitation function as obtained from the RHIC beam energy scan in Sec. 3.4; predictions 
for dielectron and dimuon spectra at current (2.76 ATeV) and future (5.5 ATeV) LHC 
energies are presented in Sec. 3.5. We end with a summary and outlook in Sec. 4. 

2. THERMAL DILEPTON RATES IN QCD MATTER 

The basic quantity for connecting calculations of the electromagnetic (EM) spectral 
function in QCD matter to measurements of dileptons in heavy-ion collisions is their 
thermal emission rate; per unit phase space it can be written as 

^-^/W^BMfM,^), (1) 

where L(M) is a lepton phase-space factor (=1 for vanishing lepton mass), f B denotes 
the thermal Bose distribution, and qo = \JM 2 + q 2 is the energy of the lepton pair (or 
virtual photon) in terms of its invariant mass and 3-momentum. As mentioned above, this 
observable is unique in its direct access to an in-medium spectral function of the formed 
system, namely in the vector (or EM) channel, LuIIem = |^M v ^ m ^EM- ^ * s defined via 
the the correlation function of the EM current, j'g M , as transported by the electric-charge 
carriers in the system. In quark basis, the EM current is given by the charge-weighted 
sum over flavor, 

Jwa = E e q qfq, (2) 

q=u,d,s 

while in hadronic basis it is in good approximation given by the vector-meson fields, 

V=p,(0,<j) & v 

known as vector-dominance model (VDM). Since the significance of thermal dilepton 
radiation is limited to masses below the J/\j/ mass, M<3 GeV, we will focus on the 
light- and strange-quark sector in this article. 

In the vacuum, the EM spectral function is well known from the e + e~ annihilation 
cross section into hadrons, usually quoted relative to the annihilation into dimuons as 
the ratio R = — ^ImTlEM, c ^- Fig- 1- It illustrates that the nonperturbative hadronic 
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FIGURE 1. Compilation of experimental data for the ratio, R, of cross sections for e + e — > hadrons 
over e + e~ — >■ jJ. + fJ.~ , as a function of invariant mass y/s = M. Figure taken from Ref. [17]. 



description in terms of VDM works well in the low-mass region (LMR), M< 1 GeV, 
while the perturbative partonic description appears to apply for M > 1.5 GeV. Thus, in 
URHICs, dilepton spectra in the LMR are ideally suited to study the properties of vector 
mesons in the medium. A central question is if and how these medium modifications can 
signal (the approach to) deconfinement and the restoration of the dynamical breaking 
of chiral symmetry (DBCS). After all, confinement and DBCS govern the properties of 
hadrons in vacuum. On the other hand, at masses M > 1 .5 GeV, the perturbative nature of 
the EM spectral function suggests that in-medium modifications are suppressed, coming 
in as corrections in powers of T/M and a s . In this case, invariant-mass spectra of thermal 
radiation become an excellent measure for the prevalent temperatures of the produced 
system, free from blue shifts due to the medium expansion which strongly affect p t 
spectra. 



2.1. Hadronic Matter 



Over the last two decades, broad efforts have been undertaken to evaluate the medium 
modifications of the p-meson. The latter dominates in the EM spectral function over the 
CO by about a factor of 10 (the appears to be rather protected from hadronic medium 
effects, presumably due to the OZI rule, at least for its coupling to baryons). Recent 
overviews of these efforts can be found, e.g., in Refs. [10, 18, 19]. Most approaches uti- 
lize effective hadronic (chiral) Lagrangians and apply them in diagrammatic many -body 
theory involving thermal (or density) loop corrections. The generic outcome is that of a 
substantial broadening of the p's spectral shape, with little mass shift (in a heat bath, chi- 
ral symmetry protects the p from mass shifts at order ff{T 2 ) [20]). The magnitude of the 
p's in-medium width (and/or its precise spectral shape) varies in different calculations, 
but the discrepancies can be mostly traced back to the differing contributions taken into 
account in the Lagrangian (e.g., the set of baryon and/or meson resonance excitations, or 
medium effects in the p's pion cloud). Similar findings arise when utilizing empirically 
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FIGURE 2. Dilepton rates from hot QCD matter. Left panel: results from the effective hadronic La- 
grangian approach for the in-medium p spectral function (solid line) at a temperature and chemical 
potential characteristic for chemical freezeout at full SPS energy; the effects of in-medium pion-cloud 
(long-dashed line), baryon resonances (dash-dotted line) and meson resonances (short-dashed line) are 
shown separately along with the rate based on the vacuum spectral function (dotted line). Right panel: 
Comparison of free and in-medium hadronic and partonic calculations in the isovector channel at temper- 
ature T = 170 MeV and small baryon chemical potential characteristic for RHIC and LHC conditions; the 
free and in-medium hadronic rates are based on Refs. [23, 24] . 



extracted on-shell p-meson scattering amplitudes off hadrons in linear-density approx- 
imation [21]. Since these calculations are restricted to resonances above the nominal 
pN (or pn) threshold, quantitative differences to many-body (or field-theoretic) ap- 
proaches may arise. In particular, the latter account for subthreshold excitations, e.g., 
p +N — >• AT* (1520), which induce additional broadening and associated enhancement 
of the low-mass part in the p spectral function (also causing marked deviations from 
a Breit-Wigner shape). Appreciable mass shifts are typically found in mean-field ap- 
proximations (due to large in-medium scalar fields), or in calculations where the bare 
parameters of the underlying Lagrangian are allowed to be temperature dependent [22]. 
An example for dilepton rates following from a p spectral function calculated in hot 
and dense hadronic matter at SPS energies is shown in the left panel of Fig. 2. The EM 
spectral function follows from the p-meson using VDM, Eq. (3), although corrections 
to VDM are necessary for quantitative descriptions of the EM couplings in the baryon 
sector [25, 26]. When extrapolated to temperatures around T pc , the resonance peak has 
essentially vanished leading to a structureless emission rate with a large enhancement in 
the mass region below the free p mass. The decomposition of the rate into in-medium 
self-energy contributions illustrates the important role of the pion cloud modifications, 
and of multiple low-energy excitations below the free p mass, e.g., resonance-hole 
BN , i.e., p +N — > B for off-shell p-mesons. The hadronic medium effects are slightly 
reduced at collider energies (right panel of Fig. 2), where a faint resonance structure 
appears to survive at around T pc . A recent calculation in a similar framework, combing 
thermal field theory with effective hadron Lagrangians [27] , and including both finite- 
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FIGURE 3. Left panel: Euclidean correlators of the EM current as computed in quenched thermal lattice 
QCD (data points) [14], compared to results from integrating hadronic spectral functions using Eq. 5 
without (dashed green line) and with in-medium effects (red lines, encompassing a free and in-medium 
continuum threshold) [31]. Right panel: Vector-isovector spectral functions at vanishing 3-momentum 
corresponding to the euclidean correlators in the left panel: in vacuum (green dashed line), in hadronic 
matter calculated from many -body theory at T= 180 MeV [13] (red solid line) and in a gluon plasma at 
1.4 T c extracted from thermal lattice-QCD (black solid line) [14]. 



temperature and -density contributions to the p selfenergy through baryon and meson 
resonances, shows fair agreement with the results shown in the left panel of Fig. 2. 



2.2. Quark-Gluon Plasma 

In a perturbative QGP (pQGP), the leading-order mechanism of dilepton production 
is EM quark-anti-quark annihilation as following from a free quark current in Eq. (2). 
The underlying EM spectral function is essentially given by the "naive quark model" 
curve in Fig. 1, extended all the way to vanishing mass, 
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where Cem = Y,q=u,d,s e l ( an additional phase-space factor occurs if finite current quark 
masses are included). Finite-temperature corrections at the same order in a s are induced 
by a thermal blocking factor which also produces a nontrivial three-momentum depen- 
dence [28] (for q=0 and H q =0 it simply amounts to a factor [1 — 2f F (qo/2)] where f F is 
the thermal Fermi distribution). The pertinent dilepton rate is structureless as shown by 
the dashed curve in the right panel of Fig. 2 (after 3-momentum integration). 

A consistent implementation of a s corrections in a thermal QGP at vanishing quark 
chemical potential has been achieved by resumming the hard-thermal-loop (HTL) ac- 
tion [29]. Quarks and gluons acquire thermal masses m'l 1 ~ gT, but Bremsstrahlungs- 
type contributions lead to a marked enhancement of the rate over the LO pQCD results, 
cf. the dash-dotted line in the right panel of Fig. 2. 



Recent progress in calculating dilepton rates nonperturbatively using thermal lattice 
QCD (1QCD) has been reported in Refs. [14, 15, 30]. The basic quantity computed 
in these simulations is the euclidean-time correlation function which is related to the 
spectral function, py = —HmU], via 

Uv ^ qJ) = .1 2^ Mq °^ T) sinh[, /2r] • (5) 

Pertinent results for Yly, normalized to its free (non-interacting) pQGP limit, have been 
obtained in quenched QCD for T=l.45T c at vanishing q (in which case M=qo), as shown 
by the data points in the left panel of Fig. 3. At small x, corresponding to large energies 
in the spectral function, this ratio tends to one as expected for the perturbative limit. For 
larger x, a significant enhancement develops which is associated with a corresponding 
enhancement in the low-energy (or low-mass) regime of the spectral function (and thus 
the dilepton rate). This enhancement can be better quantified by making an ansatz for 
the spectral function as done in Ref. [14] in terms of a low-energy Breit-Wigner part 
plus a perturbative continuum, 

P?(«o) = 5bw ^ 2 /4 + ^jp(l + K)^tanh(W4r) (6) 

(note that tanh(<5r /4r) = 1 - 2f F (q /2)). The strength (S B w) and width (T) of the Breit- 
Wigner, as well as a perturbative a s correction (k), are fit parameters. The large- x en- 
hancement in the correlator generates an appreciable low-energy enhancement in the 
spectral function, cf. right panel of Fig. 3. The zero-energy limit of the spectral func- 
tion defines a transport coefficient, the electric conductivity, Oem = gli m # ->-o(py/(7o)- 
Similar to the viscosity or heavy-quark diffusion coefficient, a small value for Oem> im- 
plied by a large value for T, indicates a strong coupling of the medium; e.g., in pQCD, 
Oem <* T /aj [32]. The results for the dilepton rate (or spectral function) at a smaller 
temperature of 1.1 T c are found to be very similar to the ones at 1.45 T c [30], suggesting 
a weak temperature dependence in this regime. Note, however, that the phase tansition 
in quenched QCD is of first order, i.e., a stronger variation could occur when going close 
to or below T c . On the other hand, recent results for two-flavor QCD [15] also indicate 
rather structureless spectral functions as found in quenched QCD. Ultimately, at suffi- 
ciently small temperatures, the lattice computations should recover a p -meson resonance 
peak; it will be interesting to see at which temperatures this occurs. 

For practical applications, a finite 3-momentum dependence of the 1QCD dilepton 
rate is needed, which is currently not available from the simulations. We here propose a 
"minimal" construction which is based on a matching to the 3-momentum dependence 
obtained from the leading-order pQCD photon rate [33]. The latter reads 

^dRr = _^ Imnr(M = ^ )/Wr) 



d 3 q K 



§S^»( 1+ ££'- 



The idea is now to adopt the transverse part of the EM spectral function as given by 
Eq. (7) for the 3-momentum dependence of the spectral function in Eq. (6) by replacing 
the Breit-Wigner part with it, i.e., 

-Imn r = ^^M 2 (l + 2K 2 a s ^F(M 2 )ln(l + ^ q ^]) . (8) 



CemN c 
12k 



: -M 2 (l + Ql AT (M iq )) (9) 



(for simplicity the small K and Pauli-blocking corrections in the annihilation term have 
been neglected). An additional formfactor, F(M 2 ) = A 2 / (A 2 +M 2 ), has been introduced 
to better reproduce the mass dependence of the 1QCD spectral function at vanishing 3- 
momentum. The depth and location of the minimum of p" /qo (as plotted in the right 
panel of Fig. 3) can be well reproduced by using A = 2T . It turns out that the electric 
conductivity following from this construction, Oem/T — 0.36Cem> is surprisingly close 
to the one estimated from the fit ansatz (6) to lattice data, (7em/T — (0.37 ±0.01 )Cem 
(also compatible with Ref. [34]; the systematic uncertainty in the lattice result, due to 
variations in the ansatz, is significantly larger). Finally, care has to be taken to include 
a finite longitudinal part which develops in the timelike regime. Here we employ a 
dependence that follows, e.g., from standard constructions of gauge-invariant 5-wave 
p-baryon interactions, yielding n^ = (M 2 /qfyllj [35]. Thus, we finally have 

Q'lat = \ (IQIat + Q L lat) = \qIat (2 + ^2) ■ (10) 

Let us compare the 1QCD results with those from hadronic calculations. In the vicinity 
of the phase transition temperature, the "melting" of the in-medium p spectral function 
suggests a smooth transition from its prominent resonance structure in vacuum to the 
rather structureless shape extracted from 1QCD, cf. right panel of Fig. 3. It would 
clearly be of interest to extract the conductivity from the hadronic calculations, which 
currently is not well resolved from the q=0, qo — > limit of the spectral function. The 
mutual approach of the nonperturbative hadronic and 1QCD spectral functions is also 
exhibited in the 3-momentum integrated dilepton rate shown in the right panel of Fig. 2, 
especially when compared the different shapes of the leading-order pQCD and vacuum 
hadronic rates. Arguably, the in-medium hadronic rate still shows an indication of a 
broad resonance. A smooth matching of the rates from above and below T pc might 
therefore require some additional medium effects in the hot and dense hadronic medium, 
and/or the emergence of resonance correlations in the qq correlator in the QGP. Unless 
otherwise noted, the thermal emission rates used in the calculations of dilepton spectra 
discussed below will be based on the in-medium hadronic rates of Ref. [24] and the 
lQCD-inspired QGP rates [14], extended to finite 3-momentum as constructed above. 



3. DILEPTON SPECTRA AT RHIC AND LHC 

To calculate dilepton mass and transverse-momentum (q t ) momentum spectra suitable 
for comparison to data in heavy-ion collisions, the thermal rates of hadronic matter and 
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FIGURE 4. Time evolution of fireball temperature in semicentral Au-Au(-y/i=0.2 GeV) collisions at 
RHIC within (the central cell of) ideal hydrodynamics (left) [40] and an expanding fireball model (right). 
The dashed green and dotted blue line in the left panel are to be compared the dashed green and solid blue 
line in the right panel, respectively. 



QGP have to be integrated over a realistic space-time evolution of the AA reaction, 
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In addition to this thermal yield, non-thermal sources have to be considered, e.g., primor- 
dial Drell-Yan annihilation and electromagnetic final-state decays of long-lived hadrons. 
We will briefly discuss space-time evolutions in Sec. 3.1 and non-thermal sources in 
Sec. 3.2, before proceeding to a more detailed discussion of thermal spectra and com- 
parisons to data, as available, in Sees. 3.3, 3.4 and 3.5 for full RHIC energy, the beam- 
energy scan and LHC, respectively. 



3.1. Medium Expansion 



The natural framework to carry out the space-time integral over the dilepton rate in 
URHICs is relativistic hydrodynamics. The application of this approach to AA collisions 
at RHIC and LHC works well to describe bulk hadron observables (e.g., p t spectra 
and elliptic flow) up to momenta of p t ~ 2 — 3 GeV, which typically comprises more 
than 90% of the total yields. Some uncertainties remain as to, e.g., the precise initial 
conditions at thermalization, viscous corrections, or the treatment of the late stages 
where the medium becomes dilute and the hadrons decouple (see, e.g., Ref. [36] for 
a recent review). Another key ingredient is the equation of state (EoS) of the medium, 
e{P), which drives its collective expansion. The left panel of Fig. 4 illustrates the effects 
of updating a previously employed bag-model EoS (a quasiparticle QGP connected to a 
hadron resonance gas via a first-order phase transition) [37 ] by a recent parametrization 
of a nonperturbative QGP EoS from 1QCD data [38, 39] (continuously matched to a 



hadron-resonance gas at 7p C =170MeV) [40]: within a 2+1-D ideal hydro calculation 
the most notable change is a significant increase of the temperature (at fixed entropy 
density) in the regime just above the transition temperature (up to ca. 30 MeV at the 
formerly onset of the 1 . order transition). Together with the fact that the hadronic portion 
of the formerly mixed phase is now entirely associated with the QGP, this will lead to 
an increase (decrease) of the QGP (hadronic) contribution to EM radiation relative to 
the first-order scenario. In addition, the harder lattice EoS induces a stronger expansion 
leading to a slightly faster cooling and thus reduction in the lifetime by about 5%. The 
latter effect becomes more pronounced if the initial conditions to the hydrodynamic 
evolution are modified, e.g., by introducing a more compact spatial profile (creating 
larger gradients) and/or an initial transverse flow (associated with interactions prior 
to the thermalization time, To) [40], cf. the solid line in the left panel of Fig. 4. The 
resulting more violent expansion plays an important role in a realistic description of the 
HBT radii of the system [41]. More importantly in the dilepton context, it can reduce 
the fireball lifetime by up to 20% or so. More simplistic parametrizations of the space- 
time evolution of AA collision have been attempted with longitudinally and transversely 
expanding fireballs. With an appropriate choice of the transverse acceleration (in all 
applications below it is taken as a r =0.12/fm at the surface), an approximate reproduction 
of the basic features (timescales and radial flow) of hydrodynamic evolutions can be 
achieved, see right panel of Fig. 4. For the most part, the results shown in the remainder 
of this article will be based on such simplified fireball parametrizations, utilizing the 
equation of state of Ref. [40] where a parametrization of 1QCD results is matched 
with a hadron resonance gas at r pc =170MeV and subsequent chemical freezeout at 
r ch =160MeV (see also Ref. [42]). We note that the use of this EoS, together with the 
lQCD-based QGP emission rates, constitutes an update of our earlier calculations [43] 
where a quasiparticle bag-model EoS was employed in connection with HTL rates in 
the QGP. We have checked that the previous level of agreement with the acceptance- 
corrected NA60 spectra is maintained, which is essentially due to the duality of the 
QGP and hadronic rates around T pc (a more detailed account in the context of the SPS 
dilepton data will be given elsewhere [45]). For our discussion of collider energies 
below, the initialization (or thermalization times) are chosen at 0.33 fm/c at full RHIC 
energy (increasing smoothly to 1 fm/c at v / 5=20GeV) and 0.2 fm/c in the LHC regime. 
This results in initial temperatures of 225 MeV and 330 MeV in minimum-bias (MB) 
Au-Au collisions at 20 and 200 GeV, respectively, increasing to ~380MeV in central 
Au-Au(200GeV), and ~560(620) MeV in central Pb-Pb at 2.76(5.5) ATeV (these values 
differ slightly from previous calculations with a quasiparticle EoS; they are also sensitive 
to the initial spatial profile, cf. left panel of Fig. 4. However, for our main objective of 
calculating low-mass dilepton spectra the initial temperature has little impact). 



3.2. Nonthermal Sources 

In addition to thermal radiation from the locally thermalized medium, dilepton emis- 
sion in URHICs can arise from interactions prior to thermalization (e.g., Drell-Yan an- 
nihilation) and from EM decays of long-lived hadrons after the fireball has decoupled 



(e.g., Dalitz decays n , r\ — > yl + l~ or ft), — > l + l~). Furthermore, paralleling the struc- 
ture in hadronic spectra, a non-thermal component from hard production will feed into 
dilepton spectra, e.g., via Bremsstrahlung from hard partons traversing the medium [44] 
or decays of both short- and long-lived hadrons which have not thermalized with the bulk 
(e.g., "hard" p-mesons or long-lived EM final-state decays). Hadronic final-state decays 
(including the double semileptonic decay of two heavy-flavor hadrons originating from 
a cc or bb pair produced together in the same hard process) are commonly referred to 
as the "cocktail", which is routinely evaluated by the experimental collaborations using 
the vacuum properties of each hadron with p t spectra based on measured spectra, or 
appropriately extrapolated using thermal blast-wave models. In URHICs, the notion of 
the cocktail becomes problematic for short-lived resonances whose lifetime is compara- 
ble to the duration of the freezeout process of the fireball (e.g., for p, A, etc.). In their 
case a better approximation is presumably to run the fireball an additional ~1 fm/c to 
treat their final-decay contribution as thermal radiation including medium effects. How- 
ever, care has to be taken in evaluating their dilepton p f -spectra, as the latter are slightly 
different for thermal radiation and final-state decays (cf. Ref. [43] for a discussion and 
implementation of this point). For light hadrons at low p t , the cocktail scales with the 
total number of charged particles, N c h, at given collision energy and centrality, while for 
hard processes a collision-number scaling <* N co n is in order (and compatible with ex- 
periment where measured, modulo the effects of "jet quenching"). The notion of "excess 
dileptons" is defined as any additional radiation observed over the cocktail, sometimes 
quantified as an "enhancement factor" in a certain invariant-mass range. The excess radi- 
ation is then most naturally associated with thermal radiation, given the usual limitations 
where hard process take over, i.e., M, q t < 2 — 3 GeV. 



3.3. RHIC-200 

We start our discussion of low-mass dilepton spectra at full RHIC energy where 
most of the current experimental information at collider energies is available, from both 
PHENIX [46] and STAR [47] measurements. 



3.3.1. Invariant-Mass Spectra 

Figure 5 shows the comparison of thermal fireball calculations with low-mass spec- 
tra from STAR [47]. As compared to earlier calculations with a bag-model EoS [13], 
the use of lQCD-EoS and emission rates for the QGP enhances the pertinent yield sig- 
nificantly. It is now comparable to the in-medium hadronic contribution for masses be- 
low M ~ 0.3 GeV, and takes over in the intermediate-mass region (M> 1.1 GeV). The 
hadronic part of the thermal yield remains prevalent in a wide range around the free p 
mass, with a broad resonance structure and appreciable contributions from Ait annihila- 
tion for M > 0.9 GeV. Upon adding the thermal yield to the final-state decay cocktail by 
the STAR collaboration (without p decay), the MB data are well described. For the cen- 
tral data, a slight overestimate around M ~ 0.2 GeV and around the ft) peak is found. A 
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FIGURE 5. Dilepton invariant-mass spectra in Au-Au(200 AGeV) for minimum-bias (left) and central 
(right) collisions. Theoretical calculations for thermal radiation from a nonperturbative QGP and in- 
medium hadronic spectral functions are compared to STAR data [47, 48]. 



similar description [49] of the STAR data arises in a viscous hydrodynamic description 
of the medium using the p spectral function from on-shell scattering amplitudes [21] 
(see also Ref. [50]), and in the parton-hadron string dynamics transport approach using 
a schematic temperature- and density-dependent broadening in a Breit-Wigner approxi- 
mation of the p spectral function [51]. More studies are needed to discern the sensitivity 
of the data to the in-medium spectral shape, as the latter significantly varies in the differ- 
ent approaches. For the PHENIX data (not shown), the enhancement as recently reported 
in Ref. [52] for non-central collisions (carried out with the hadron-blind detector (HBD) 
upgrade) agrees with earlier measurements [46], is consistent with the STAR data, and 
thus should also agree with theory. For the most central Au-Au data, however, a large 
enhancement was reported in Ref. [46] which is well above theoretical calculations with 
broad spectral functions [13, 53, 51, 54], even in the MB sample. More "exotic" expla- 
nations of this effect, which did not figure at the SPS, e.g., within a Bose-condensed 
like glasma in the pre-equilibrium stages [55], have been put forward to explain the 
"PHENIX puzzle". However, it is essential to first resolve the discrepancy on the exper- 
imental side, which is anticipated with the HBD measurement for central collisions. 

To quantify the centrality dependence of the thermal radiation (or excess) yield, one 
commonly introduces an exponent a c as YM(N c h)/N c h = C/V c ", which describes how the 
excess (or thermal) yield at a given mass (or mass range) scales relative to the charged- 
particle multiplicity. For full RHIC energy, the theoretical calculation gives a c ~ 0.45 
(with a ca. 10% error), similar to what had been found for integrated thermal photon 
yields in Ref . [56]. 
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3.3.2. Transverse-Momentum Dependencies 

When corrected for acceptance, invariant-mass spectra are unaffected by any blue- 
shift of the expanding medium, which renders them a pristine probe for in-medium 
spectral modifications. However, the different collective flow associated with different 
sources may be helpful in discriminating them by investigating their q t spectra, see, e.g., 
Refs. [57, 58, 28, 59, 60, 61, 62]. As is well-known from the observed final-state hadron 
spectra, particles of larger mass experience a larger blue- shift due to collective motion 
in the expanding medium than lighter particles. Schematically, this can be represented 
by an effective slope parameter which for sufficiently large masses takes an approximate 
form of r e ff = T +M/3 2 where T and /3 denote the temperature and average expansion 
velocity of the emitting source. Dileptons are well suited to systematically scan the mass 
dependence of T e ff by studying q t spectra for different mass bins (provided the data have 
sufficient statistics). At the SPS this has been done by the NA60 collaboration [63], 
who found a gradual increase in the slope from the dimuon threshold to the p mass 
characteristic for a source of hadronic origin (aka in-medium p mesons), a maximum 
around the p mass (late p decays), followed by a decrease and leveling off in the inter- 
mediate mass region (IMR, M > 1 GeV) indicative for early emission at temperatures 
T ~ 170 - 200 MeV (where at the SPS the collective flow is still small). 

Figure 6 shows the q t spectra for thermal radiation from hadronic matter and QGP in 
MB Au-Au(200 AGeV) in two typical mass regions where either of the two sources 
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FIGURE 7. Dilepton invariant mass-spectra in different bins of transverse momentum from thermal 
radiation of QGP (dash-doted line) and hadronic matter (solid line: in-medium, dashed line: vacuum 
spectral function) in MB Au-Au(200 AGeV) collisions; experimental acceptance as in Figs. 5 and 6. 



dominates. In the low-mass region (LMR), both sources have a surprisingly similar 
slope (T s i p e ~ 280 MeV), reiterating the notion that the emission is from mostly around 
T c where the slope of both sources is comparable (also recall from Fig. 5 that in the 
mass window M=0. 3-0.7 GeV the QGP emission is largest at the lower mass end, while 
the hadronic one is more weighted toward the higher end). For definiteness, assuming 
r=170MeV and M=0.5 GeV, one finds /3 ~ 0.45 - 0.5, which is right in the expected 
range [40]. On the other hand, in the IMR, where the QGP dominates, the hadronic slope 
has significantly increased to ca. 360 MeV due to the larger mass in the collective-flow 
term. On the other hand, the slope of the QGP emission has only slightly increased over 
the LMR, indicating that the increase in mass in the flow-term is essentially offset by an 
earlier average emission temperature, as expected at higher mass (for hadronic emission, 
the maximal temperature is obviously limited by T pc ). Taken together, this implies that at 
RHIC the effective slope of the total thermal radiation in the IMR significantly exceeds 
the one in the LMR, contrary to what has been observed at SPS. Together with blue-shift 
free temperature measurements from slopes in invariant-mass spectra, this provides a 
powerful tool for disentangling collective and thermal properties through EM radiation 
from the medium. 
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FIGURE 8. Left panel: Inclusive elliptic flow of thermal dileptons in 0-20% central Au-Au(200 AGeV) 
collisions, calculated within an ideal hydrodynamic model with lattice EoS using lQCD-based QGP and 
medium-modified hadronic rates [67]. Right panel: Dielectron-V2 measured by STAR in MB Au-Au [68], 
including the cocktail contribution; the latter has been simulated by STAR and is shown separately by the 
solid histogram. 



Alternatively, one can investigate the mass spectra in different momentum bins, pos- 
sibly revealing a q t -dependence of the spetral shape, as was done for both e + e~ data 
in Pb-Au [64] and jl + jl in In-In [63] at SPS. Calculations for thermal radiation in 
Au-Au at full RHIC energy are shown in Fig. 7 for four bins from q t =0-2 GeV. One 
indeed recognizes that the p resonance structure becomes more pronounced as trans- 
verse momentum is increased. In the lowest bin, there is no low-mass crossing with the 
QGP contribution anymore; the minimum structure around M ~ 0.2 GeV is caused by 
the experimental acceptance, specifically the single electron p e t > 0.2 GeV, which for 
vanishing q t suppresses all dilepton yields below M ~ 2p e t ,min = 0.4 GeV. 



3.3.3. Elliptic Flow 

Another promising observable to diagnose the collectivity, and thus the origin of the 
EM emission source, is its elliptic flow [65, 62, 66]. The latter is particularly useful to 
discriminate early from late(r) thermal emission sources: contrary to the slope parameter, 
which is subject to an interplay of decreasing temperature and increasing flow, the 
medium's ellipticity is genuinely small (large) in the early (later) phases. The left panel 
of Fig. 8 shows hydrodynamic calculations of the inclusive thermal dilepton V2 as a 
function of invariant mass (using the same emission rates and EoS as in previous figures 
shown here) [67]. One nicely recognizes a broad maximum structure around the p mass, 
indicative for predominantly later emission source in the vicinity of its vacuum mass, a 
characteristic mass dependence (together with an increasing QGP fraction) below, and 
a transition to a dominant QGP fraction with reduced V2 above. All these features are 
essentially paralleling the mass dependence of the slope parameter at SPS, while the 
latter exhibits a marked increase at RHIC in the IMR due to the increased radial flow in 
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FIGURE 9. Low-mass dilepton spectra as measured by STAR in the RHIC beam-energy scan [68]; MB 
spectra are compared to theoretical predictions for the in- medium hadronic + QGP radiation, added to the 
cocktail contribution. 



the QGP and early hadronic phase. Rather similar results are obtained in hydrodynamic 
calculations with in-medium spectral functions from on-shell scattering amplitudes [49]. 
When using a less pronounced in-medium broadening, the peak structure in V2(M) 
tends to become narrower [65, 62, 66]. First measurements of the dilepton-V2 have been 
presented by STAR [68], see right panel of Fig. 8. The shape of the data is not unlike 
in the theoretical calculations, while it is also consistent with the simulated cocktail 
contribution. Note that the total V2 is essentially a weighted sum of cocktail and excess 
radiation. Thus, if the total V2 were to agree with the cocktail, it implies that the V2 of the 
excess radiation is as large as that of the cocktail. Clearly, future V2 measurements with 
improved accuracy will be a rich source of information. 

Significant V2 measurements of EM excess radiation have recently been reported in 
the M=0 limit, i.e., for direct photons, by both PHENIX [69, 70] and ALICE [71, 72]. 
A rather large vj signal has been observed in both experiments [70, 72], suggestive for 
rather late emission [73] (see also Refs. [74, 75, 76, 77]). In addition, the effective slope 
parameters of the excess radiation have been extracted, T e ff = 219±27MeV [69] at 
RHIC-200 and 304 ± 51 MeV at LHC-2760 [71], which are rather soft once blue-shift 
effects are accounted for. In fact, these slopes may not be inconsistent with the trends in 
the LMR dileptons when going from RHIC (Fig. 6 above) to LHC (Fig. 12 below). This 
would corroborate their main origin from around r pc . 



3.4. RHIC Beam Energy Scan 



A central question for studying QCD phase structure is how the spectral properties 
of excitations behave as a function of chemical potential and temperature. With the 
EM (or vector) spectral function being the only one directly accessible via dileptons, 
systematic measurements as a function of beam energy are mandatory. At fixed target 
energies, this is being addressed by the current and future HADES efforts (£'beam=l- 



10 AGeV) [78, 79], by CBM for £ bea m(Au) up to ~35 AGeV [2], and has been measured 
at SPS energies at £ beam =158 AGeV, as well as in a CERES run at 40 AGeV [80]. 

At collider energies, a first systematic study of the excitation function of dilepton 
spectra has been conducted by STAR [68] as part of the beam-energy scan program 
at RHIC. The low-mass excess radiation develops rather smoothly when going down 
form y/SNN=200 GeV via 62GeV to 20GeV, cf. Fig. 9. Closer inspection reveals that 
the enhancement factor of excess radiation over cocktail in the region below the p mass 
increases as the energy is reduced [68]. An indication of a similar trend was observed 
when comparing the CERES measurements in Pb-Au at v /i7v^v=17.3 GeV and 8.8 GeV. 
Theoretically, this can be understood by the importance of baryons in the generation of 
medium effects [26], specifically the low-mass enhancement in the p spectral function. 
These medium effects become stronger as the beam energy is reduced since the hadronic 
medium close to T pc becomes increasingly baryon rich. At the same time, the cocktail 
contributions, which are mostly made up by meson decays, decrease. The hadronic in- 
medium effects are expected to play a key role in the dilepton excess even at collider 
energies. The comparison with the STAR excitation function supports the interpretation 
of the excess radiation as originating from a melting p resonance in the vicinity of T pc . 

A major objective of the beam-energy program is the search for a critical point. 
One of the main effects associated with a second-order endpoint is the critical slowing 
down of relaxation rates due to the increase in the correlation length in the system. 
For the medium expansion in URHICs, this may imply an "anomalous" increase in 
the lifetime of the interacting fireball. If this is so, dileptons may be an ideal tool 
to detect this phenomenon, since their total yield (as quantified by their enhancement 
factor), is directly proportional to the duration of emission. The NA60 data have shown 
that a pertinent lifetime measurement can be carried out with an uncertainty of about 
±1 fm/c [43]. In the calculations shown in Fig. 9 no critical slowing down has been 
assumed; as a result, the average lifetime in MB Au-Au collisions increases smoothly 
from ca. 8 to 10 fm/c. Thus, if a critical point were to exist and lead to a, say, 30% 
increase in the lifetime in a reasonably localized range of beam energies, dilepton yields 
ought to be able to detect this signature. This signal would further benefit from the fact 
that the prevalent radiation arises from around T pc where the largest effect from the 
slowing down is expected. 



3.5. LHC 

The previous section raises the question whether the smooth excitation function of 
dilepton invariant-mass spectra in the RHIC regime will continue to LHC energies, 
which increase by another factor of ~20. On the other hand, the dilepton q t spectra, 
especially their inverse-slope parameters, indicate an appreciable variation from SPS to 
RHIC, increasing from ca. 220 to 280 MeV in the LMR, and, more pronounced, from 
ca. 210 to about 320 MeV in the IMR. This is a direct consequence of the much stronger 
(longer) development of collective flow in the QGP phase of the fireball evolution. 
This trend will continue at the LHC, as we will see below. In the following section 
(3.5.1), we will first discuss the dielectron channel at LHC, and highlight the excellent 
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FIGURE 10. Dielectron invariant-mass spectra from thermal radiation in 0-40% central Pb- 
Pb(2.76ATeV) (left panel) and 0-10% central Pb-Pb(5.5ATeV) (right panel), including single-electron 
cuts to simulate the ALICE acceptance. Hadronic (with in-medium or vacuum EM spectral function) and 
QGP contributions are shown separately along with the sum of in-medium hadronic plus QGP. Here and 
in the following LHC plots, both vacuum and in-medium hadronic emission rates in the LMR have been 
supplemented with the vacuum spectral function in the LMR, i.e., no in-medium effects due to chiral 
mixing have been included (for all RHIC calculations shown in the previous sections full chiral mixing 
was included). 



experimental capabilities that are anticipated with a planned major upgrade program of 
the ALICE detector [81]. In addition, ALICE can measure in the dimuon channel, albeit 
with somewhat more restrictive cuts whose impact will be illustrated in Sec. 3.5.2. 



3.5.1. Dielectrons 



The invariant-mass spectra of thermal radiation at LHC energies show a very similar 
shape and hadronic/QGP composition as at RHIC energy, see Fig. 10. This may not 
be surprising given the virtually identical in-medium hadronic and QGP rates along the 
thermodynamic trajectories at RHIC and LHC (where /ig <C T at chemical freezeout). It 
also implies that the thermal radiation into the LMR is still dominated by temperatures 
around T pc , with little (if any) sensitivity to the earliest phases. The total yield, on the 
other hand, increases substantially due to the much larger fireball volumes created by 
the larger multiplicities. More quantitatively, the (A^h-normalized) enhancement around, 
e.g., M=0.4 GeV, approximately scales as A^ £ with Ofe ~ 0.8 relative to central Au-Au at 
full RHIC energy. This is a significantly stronger increase than the centrality dependent 
enhancement at fixed collision energy, a c ~ 0.45 as quoted in Sec. 3.3.1. 

Detailed simulation studies of a proposed major upgrade of the ALICE detector have 
been conducted in the context of a pertinent letter of intent [81]. The final results after 
subtraction of uncorrected (combinatorial) background are summarized in Fig. 11, 
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FIGURE 11. Simulations of dielectron invariant-mass spectra in Pb-Pb(5.5 ATeV) collisions assuming 
the thermal spectra shown in the left panel of Fig. 10 as the excess signal [81, 82]. In addition to the 
accpetance cuts on single-electron rapidity and momentum, pair efficiency and displaced vertex cuts 
are included here. Left panel: Invariant-mass spectra after subtraction of combinatorial background; in 
addition to the thermal signal the simulated data contain the hadronic cocktail and correlated open- 
charm decays. Right panel: Simulated excess spectra after subtraction of the cocktail and the open-charm 
contribution using displaced vertex cuts. 



based on an excess signal given by the thermal contributions in Fig. 10. 1 The left 
panel shows that the thermal signal is dominant for the most part of the LMR (from 
ca. 0.2-1. GeV), while in the IMR it is outshined by correlated heavy-flavor decays. 
However, the latter can be effectively dealt with using displaced vertex cuts, while 
the excellent mass resolution, combined with measured and/or inferred knowledge of 
the Dalitz spectra of n° (from charged pions), r\ (from charged kaons) and CO (from 
the direct dilepton decay), facilitate a reliable subtraction of the cocktail. The resulting 
excess spectra shown in the right panel are of a quality comparable to the NA60 data. 
This will allow for quantitative studies of the in-medium EM spectral function in the 
LMR which are critical for being able to evaluate signatures of chiral restoration (as 
discussed elsewhere, see, e.g., Refs. [83, 16]). In addition, the yield and spectral slope of 
the dominantly QGP emission in the IMR will open a pristine window on QGP lifetime 
and temperature (recall that the M spectra, which are little affected by the acceptance 
cuts in the IMR, are unaffected by blue shifts). 

Let us turn to the dilepton q t spectra at full LHC energy, displayed again for two 
mass bins representing the LMR and IMR in Fig. 12. Compared to RHIC, the LHC 



1 The thermal yields provided for the simulations were later found to contain a coding error in the author's 
implementation of the experimental acceptance; the error turns out to be rather inconsequential for the 
shape and relative composition of the signal, as a close comparison of the right panels in Fig. 10 and 1 1 
reveals; the differentail yields differ by up to 20-30%. 







_ , , , , 1 , , , , 1 , , 


1 1 ' ' ' ' 1 ' ' ' ' - 


! , 






— hadronic (0.3-0.7) - 




w 5 


K._"] ^-^445 


QGP (0.3-0.7) i 


o 




- ' ^* .^-V 


— hadronic (1.1-1.5) : 




ID 6 


. 450 '^Sv 


QGP (1.1-1.5) J 


>» 




"""'•-. ''On. 


z 


^ 




" ' *> . ^465 ' V\ 


- 


-C 




_ ~ • "vV- 




i 


io- 7 




VS. E 

VS Z 


^ 




"*" N_ N _ - 


*■ ^"\- 


w 




" s ' . 


' ^ . ^^A* 




10 8 


_ ^ 




>> 




Central Pb-Pb (5.5TeV) 




5 


io 9 


<N ,>=1930 

— ch 




1 




p t c >0.2GeV, ly c l<0.84 


* 




, „-io 


, , I , 


, 1 



2 3 

q t [GeV] 



FIGURE 12. Same as Fig. 6, but for central Pb-Pb(5.5 ATeV). 



fireball is characterized by a marked increase in QGP lifetime and associated build-up 
of transverse flow by the time the system has cooled down to T pc . Consequently, the q t 
spectra exhibit an appreciable increase in their inverse- slope parameters, by about 60% 
in the LMR (for both hadronic and QGP parts) and for the QGP part in the IMR, and 
up to 80% for the hadronic IMR radiation (recall that in a scenario with chiral mixing, 
the hadronic radiation for Af=l. 1-1.5 GeV is expected to increase by about a factor of 
2, so that its much larger slope compared to the QGP will become more significant for 
the total). 



3.5.2. Dimuons 



Low-mass dilepton measurements are also possible with ALICE in the dimuon chan- 
nel at forward rapidities, 2.5 < y^ < 4, albeit with somewhat more restrictive momentum 
cuts [84]. The charged-particle multiplicity in this rapidity range is reduced by about 
30% compared to midrapidity [85], but, at 2.76 ATeV, is still ca. 30% above central 
rapidities in central Au-Au at RHIC. 

Figure 13 illustrates the expected thermal mass spectra in central Pb-Pb(2.76 ATeV). 
For "conservative" cuts on the di-/muons (q^ > 2 GeV, pf > 0.7 GeV), their yield is 
substantially suppressed (see left panel), by about one order of magnitude, compared 
to a typical single-e cut of p e t > 0.2 GeV In addition, the spectral broadening of the 
in-medium p meson is less pronounced, a trend that was also observed in the ^-sliced 
NA60 dimuon spectra. Here it is mostly due to the suppression of medium effects at 
larger p-meson momentum relative to the heat bath, caused by hadronic formfactors 
(analogous to RHIC, recall Fig. 7). It is, in fact, mostly the pair cut which is responsible 
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FIGURE 13. Calculations of thermal dimuon invariant-mass spectra in central Pb-Pb(2.76 ATeV) col- 
lisions at forward rapidity, y=2.5-4. In the left panel, in-medium hadronic, vacuum hadronic, QGP and 
the sum of in-medium hadronic plus QGP are shown with "strong" cuts on single- and dimuon transverse 
momenta. The right panel illustrates how the total yield increases when the two cuts are relaxed. 



for the suppression, since q^ ,cu is significantly larger than 2p[ t ' cu . If the former can be 
lowered to, say, 1 GeV, the thermal yield of accepted pairs increases by about a factor 
3 in the IMR and 2 in the LMR (see dashed line in the right panel of Fig. 13). The 
LMR acceptance is now mainly limited by the single-/! cut, as the latter suppresses low- 
mass pairs whose pair momentum is not at least 2p^ ,cu (the same effect is responsible 
for the rather sharp decrease in acceptance for low-momentum electron pairs below 
M ~ 0.4 GeV in the upper left panel of Fig. 7, leading to a dip toward lower mass in the 
thermal spectra, even though the thermal rate increases approximately exponentially). 
This could be much improved by lowering the single-/! cut to, e.g., 0.5 GeV, which 
would increase the low-mass yield by about a factor of 3. At the same time, the spectral 
broadening of the p becomes more pronounced in the accepted yields, i.e., the data 
would be more sensitive to medium effects. 

4. SUMMARY AND OUTLOOK 



In this article we have given an overview of medium modifications of the electromag- 
netic spectral function under conditions expected at collider energies (high temperature 
and small baryon chemical potential), and how these medium effects manifest them- 
selves in experimental dilepton spectra at RHIC and LHC. For the emission rates from 
the hadronic phase we have focused on predictions from effective hadronic Lagrangians 
evaluated with standard many-body (or thermal field-theory) techniques; no in-medium 
changes of the parameters in the Lagrangian (masses and couplings) have been assumed. 
Since this approach turned out to describe dilepton data at the SPS well, providing 
testable predictions for upcoming measurements at RHIC and LHC is in order. As colli- 
sion energy increases the QGP occupies an increasing space-time volume of the fireball 



evolution. To improve the description of the pertinent dilepton radiation, information 
from lattice-QCD has been implemented on (i) the equation of state around and above 
T pc , and (ii) nonperturbative dilepton emission rates in the QGP. The latter have been 
"minimally" extended to finite 3-momentum to facilitate their use in calculations of ob- 
servables. Since these rates are rather similar to previously employed perturbative (HTL) 
rates, an approximate degeneracy of the in-medium hadronic and the 1QCD rates close 
to T pc still holds. This is welcome in view of the smooth cross-over transition as deduced 
from bulk properties and order parameters at jl q =0. 

The main features of the calculated thermal spectra at RHIC and LHC are as follows. 
The cross-over in the 1QCD EoS produces a noticeable increase of the QGP fraction 
of the yields (compared to a bag-model EoS), while the hadronic portion decreases (its 
former mixed-phase contribution has been swallowed by the QGP). However, due to the 
near-degeneracy of the QGP and hadronic emission rates near T pc , both the total yield 
and its spectral shape change little; the hadronic part remains prevalent in an extended 
region around the p mass at all collision energies. The very fact that an appreciable 
reshuffling of hadronic and QGP contributions from the transition region occurs, indi- 
cates that the latter is a dominant source of low-mass dileptons at both RHIC and LHC. 
This creates a favorable set-up for in-depth studies of the chiral restoration transition 
in a regime of the phase diagram where quantitative support from 1QCD computations 
for order parameters and the EM correlator is available. Current ideas of how to render 
these connections rigorous have been reported elsewhere. Phenomenologically, it turns 
out that current RHIC data for LMR dilepton spectra are consistent with the melting 
p scenario 2 , including a recent first measurement of an excitation function all the way 
down to SPS energies. However, the accuracy of the current data does not yet suffice 
to discriminate in-medium spectral functions which differ considerably in detail. These 
"details" will have to be ironed out to enable definite tests of chiral restoration through 
the EM spectral function. 

While the low-mass shape of the spectra is expected to be remarkably stable with 
collision energy, large variations are predicted in the excitation function of other dilepton 
observables. First, the total yields increase substantially with collision energy. In the 
LMR, the dependence on charged-particle rapidity density, N°^, is estimated to scale 
as OC ~ 1.8 from RHIC to LHC, significantly stronger than as function of centrality 
at fixed y^. This is, of course, a direct consequence of the longer time it takes for the 
fireball to cool down to thermal freezeout. For the RHIC beam-energy scan program, this 
opens an exciting possibility to search for a non-monotonous behavior in the fireball's 
lifetime due to a critical slowing down of the system's expansion. If the LMR radiation 
indeed emanates largely from the transition region, a slowed expansion around T c would 
take maximal advantage of this, thus rendering an "anomalous dilepton enhancement" a 
promising signature of a critical point. 

Second, the transverse-momentum spectra of thermal dileptons are expected to be- 
come much harder with collision energy, reflecting the increase in the collective ex- 
pansion generated by the QGP prior to the transition region. This "QGP barometer" 
provides a higher sensitivity than final- state hadron spectra which include the full col- 



2 With the caveat of the central Au-Au PHENIX data. 



lectivity of the hadronic evolution. The inverse-slope parameters for q t spectra in the 
LMR are expected to increase from ~220 MeV at SPS, to ~280 MeV at RHIC-200, and 
up to ~450 MeV at LHC-5500. Even larger values are reached in the IMR, although the 
situation is a bit more involved here, since (a) the QGP emission is increasingly emitted 
from earlier phases, and (b) the hadronic emission, while picking up the full effect of 
additional collectivity at T pc , becomes subleading relative to the QGR The trend in the 
LMR seems to line up with the recent slope measurements in photon excess spectra at 
RHIC and LHC. A similar connection exists for the elliptic flow; pertinent data will be 
of great interest. Invariant-mass spectra in the IMR remain the most promising observ- 
able to measure early QGP temperatures, once the the correlated heavy-flavor decays 
can be either subtracted or reliably evaluated theoretically. 

The versatility of dileptons at collider energies comprises a broad range of topics, 
ranging from chiral restoration to direct-temperature measurements, QGP collectivity, 
and fireball lifetime. Experimental efforts are well underway to exploit these, while 
sustained theoretical efforts will be required to provide thorough interpretations. 
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